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We perform density functional calculations of the equilibrium ^^Fe/^^Fe ratios for ferrous iron 
dissolved in periclase and MgSiOs perovskite at the pressures and temperatures of the Earth's 
mantle. Pressure increases the partitioning of ^^Fe into both phases by a factor of three from the 
Earth's surface to the core-mantle boundary. In ferropericlase, a large contribution to this increase 
comes from the electronic transition from high-spin to low-spin iron. Our calculations demonstrate 
that pressure-induced fractionation can play a major role in determining planetary iron-isotope 
composition. 
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The distribution of iron-isotopes between and within 
planetary bodies provides constraints on their histories 
of accretion and differentiation [U [2l [3]. A significant 
fraction of the iron in the silicate Earth is dissolved 
in periclase (Fca^Mgi-a^O) (fep) and silicate perovskite 
(FeajMgi_a^)Si03 (fepv) which make up the bulk of the 
Earth's lower mantle [4 . Uptake by fep and fepv has al- 
most certainly influenced the distribution of iron-isotopes 
in the Earth, but little is known definitively about the 
equilibrium iron isotope compositions of these phases at 
lower-mantle pressures and temperatures. Empirical es- 
timates, based on the vibrational density of states for the 
iron-sublattice obtained from Mossbauer spectroscopy 
and inelastic nuclear resonant X-ray scattering (INRXS) 
techniques, are starting to address this issue [5l[6]. 

Equilibrium isotopic signatures should become heav- 
ier with increasing pressure as bonds are compressed and 
vibrational frequencies increase. Because the ionic ra- 
dius of low-spin iron is smaller than high-spin iron, iron- 
isotope signatures also may be made heavier by the spin 
transition for Fe^+ [ZIE], which, in fep, is predicted to 
occur within the Earth's mantle [9, 10 . The reduced par- 
tition function ratio (/3), can be used to make predictions 
of the effect of pressure and electronic state on the possi- 
ble isotopic fractionations associated with fep and fepv. 
In the harmonic approximation (3 is defined by pT] : 
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where ^^(/i,/) . = hc27ruj(^h^i^ -/kT, h and / refer to the heavy 
isotope and light isotope, respectively (here ^^Fe and 
^^Fe), T is the temperature, and the product runs over all 
frequencies ooi. The equilibrium constant for isotopic ex- 
change a between phases i and j is given as aij = Pi/Pj] 
aij > 1 implies that the heavy isotope is concentrated in 
phase i. Here we use density functional electronic struc- 
ture methods to calculate p for Fe^+ in fepv and fep at 
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pressures and temperatures near that of the core-mantle 
boundary. We have found that p increases dramatically 
with pressure and is strongly affected by the high-spin to 
low-spin transition in fep. 

Density functional calculations have made important 
contributions to understanding the properties of lower 
mantle phases [HI [TOl [12]. They also have been used 
to estimate isotope fractionation factors in iron-bearing 
phases such as Fe ^' ^(aq) via harmonic vibrational fre- 
quencies [T3","T4 . Here, we use these methods on clusters 
representative of the coordination environment of Fe^+ 
in fep and fepv. The clusters can be treated with stan- 
dard all-electron quantum chemistry methods, avoiding 
potential problems of the iron pseudopotential and mak- 
ing a stronger connection with existing studies of vibra- 
tional frequencies and the spin-crossover transition in the 
quantum chemistry literature, which can guide the choice 
of exchange-correlation functionals and basis set [15]. 
Such methods allow affordable use of hybrid exchange- 
correlation functionals which can make accurate esti- 
mates of the spin-crossover energies in extensive studies 
on molecular systems [16 and have provided the best 
estimates of isotope fractionation factors [17]. The hy- 
brid functionals complement calculations using the local- 
density and generalized-gradient approximations [HI [T2] 
and studies using the Hubbard model [TO] . 

Iron atoms are assumed to vibrate in coordination en- 
vironments characteristic of fep and fepv, similar to 
the site-specific vibrational density of states approach 
adopted in Ref. [6] . The models consist of a free "core" 
surrounded by a rigid shell of oxygen atoms fixed in their 
ideal lattice positions. The vibrational modes of the core 
within the cavity that displace the iron center are as- 
sumed to approximate the vibrational modes that gov- 
ern equilibrium isotopic fractionation in the crystal. The 
cores of the clusters are shown in Figure 1. The fepv 
site consists of a single iron atom and eight coordinating 
oxide atoms. All silicon and magnesium atoms attached 
to these eight oxygens are also included in the core clus- 
ter. The core is surrounded by a layer of oxygen atoms 
completing the coordination shells of the silicon and mag- 




TABLE I: Calculated bond lengths, in A, with multiplicities 
given in parentheses, if different from unity. 



FIG. 1: Core atoms of molecular clusters used to model fe- 
periclase (fep, above) and fe-perovskite (fepv, below). The 
iron atom (brown) is centered within its coordination environ- 
ment of oxygen (red) , magnesium (yellow) and silicon (blue) . 
These core atoms are fully relaxed within a shell of oxygen 
atoms, fixed in their ideal lattice positions. This rigid outer 
shell of oxygen atoms is capped by hydrogen atoms having 
fractional charges as described in the text. 



nesium atoms. The other bonds into the outer layer of 
oxygen atoms coming from "external" magnesium and 
silicon atoms, are represented by link atoms having a 
charge equal to the Pauling bond strength (the charge di- 
vided by the coordination number) contributed into the 
oxygen atom. The link atoms are placed along the broken 
bonds one angstrom away from the oxygen atom. This 
procedure ensures a neutral, autocompensated cluster. 
The fep cluster, centered on the octahedral iron site, is 
constructed in the same manner. This simple technique 
accurately reproduces calculations of isotope fractiona- 
tion in carbonates and other oxide minerals with periodic 
boundary conditions [171 EH] • 

The electronic structure calculations, carried out with 
the PQS quantum chemistry package of Pulay, Baker, 
and co-workers (http://www.pqs-chem.com|). The rigid 
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2.478 (2) 2.212 (2) 



2.213 (2) 1.978 (2) 2.213 (6) 1.959 (6) 
2.276 (4) 2.008 (4) 



outer shell uses the 3-2 IG basis on both the oxygen atoms 
and the link atoms (which use hydrogen basis functions) . 
The optimized atoms use the 6-31 G* basis for O, Mg, 
and Si, and the m6-31G* basis for Fe [19 . The B3LYP 
exchange-correlation functional, as implemented in the 
PQS code, was used for the exchange-correlation func- 
tional. For fepv outer-shell oxygen positions were taken 
from Ref. [20] at GPa and from electronic structure 
calculations [2T| at 120 GPa. Lattice parameters for fep 
were obtained from a recent universal P-V-T equation of 
state for periclase [22 . The high-pressure lattice param- 
eters were taken at 120 GPa and at the highest available 
temperature of 3000 K, even though we evaluate Eq. 1 at 
4000 K. The low-pressure lattice parameters were taken 
at GPa and 300 K. 

For the fepv cluster, we enforce mirror symmetry 
through the center of the cluster perpendicular to the 
c crystallographic axes, giving the cluster point symme- 
try Cs' For the low-spin fep cluster we retain the T/j, 
symmetry group; this is relaxed in the high-spin cluster 
which was run with no symmetry constraints to acco- 
modate possible Jahn- Teller distortions. The positions 
of the central iron atom, its coordinating oxygen atoms, 
and the magnesium and silicon atoms attached to the co- 
ordinating oxygen atoms are optimized within the rigid 
outer shell of the cluster (oxygen atoms and link atoms) . 
Relaxed bond lengths in the primary coordination sphere 
of the iron for all clusters are given in Table 1. 

The Hessian matrix for the free atoms was calculated 
analytically by projecting out the degrees of freedom of 
the fixed atoms. Harmonic vibrational frequencies were 
calculated for the core atoms with both ^^Fe and ^^Fe, 
and used to estimate j3 in the harmonic approximation 
according to Eq. 1. The j3 defined by Eq. 1 was divided 
by the high temperature limiting value Pht according to 
standard convention [2T. 

Calculations on low-spin Fe^+ in the fepv cage in Cs 
symmetry gave a vibrational mode with an imaginary 
frequency perpendicular to the mirror plane, low- spin 
Fe^+ would therefore "rattle" unfavorably in a cage with 
Cs symmetry, qualitatively consistent with the idea the 
Fe^+ is not likely to undergo a pressure-induced high- 
spin to low-spin transition in fepv in the Earth's mantle 
[24] . We did not consider low-spin Fe^+ further in the 
fepv cluster. Fe^+ in the octahedral Si^^ site may un- 



TABLE II: Reduced partition function ratios 10^\n{/3/(3HT), 
calculated for ^^Fe amd ^'^Fe, according to Eq. 1 and as de- 
scribed in the text. 





fepv 

high-spin 

120 


fep 

high-spin low-spin 


P (GPa) 


120 120 


300 K 
4000 K" 


9.86 24.2 
0.06 0.15 


7.04 (7.8)" 16.5 8.1 24.1 (20.1)"^ 
0.04 0.10 0.05 0.15 



^taken from Ref. [6] 
^P = 109 GPa 

^the periclase lattice parameters are taken at 120 GPa and 3000 K 
from Ref. [22j but the /3 was evaluated at 4000 K in Eq. 1 



dergo spin crossover within the mantle [T?, but this is 
not considered here. 

The results of the calculations are shown in Table 2, 
given as 10Hn{P/PHT)^ for temperatures of 300 K and 
4000 K, the latter temperature being close to that of the 
core-mantle boundary. First, pressures at the core man- 
tle boundary increase the /3 factors for both phases by a 
factor of 2.5-3 relative to the Earth's surface. For fep at 
pressures near the core-mantle boundary, half of the en- 
richment can be attributed to the high-spin to low-spin 
transition in fep A much smaller difference is predicted 
at ambient pressures. Although we do not claim to be 
able to quantitatively model the energetics of the spin 
transition here, we remark that the electronic energy for 
the cluster with low-spin iron fep is approximately 1.5 eV 
above that for the cluster with high-spin iron at GPa. 
In the high-pressure cluster, however, the electronic en- 
ergy of the low-spin cluster is about 0.4 eV below that for 
high-spin cluster. The /3 values for fep given in Table 2 
are in reasonable agreement with those reported in [6 . 
The high-pressure value of P is is somewhat lower than 
our calculations. Part of the reason for this is that the 
high-pressure value of P based on INRXS measurements 
is given at 109 GPa whereas the calculations were carried 
out at ^ 120 GPa. 

The contributions to /3 for the high-spin and low-spin 
fep clusters are shown in Fig. 2. Gontributions to /3 
are spread out among seven main fractionating modes. 
Gharacteristic motions associated with these modes are 
given in Fig. 3. For both the high-spin and low-spin 
configurations, the largest contribution to /3 comes from 
the fundamental vibration of the central iron against the 
cage of surrounding oxygens and magnesium ions at 308 
(302) cm~-^ for ^^Fe(^^Fe) in low-spin fep. This is an 
accoustic mode with the coordinating oxygen atoms and 
magnesium atoms moving with the iron, but with less 
displacement. Although this mode makes the biggest 
contribution to /?, the difference between the high-spin 
and low-spin contributions is not very large. The next 
most important fractionating mode is at 491.3 (489.80) 
cm~^, an optical mode involving the iron atom vibrat- 
ing against its own foxygen coordination shell and their 
attached magnesium atoms. In the high-spin fep clus- 
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FIG. 2: Gumulative value of {(3-1) for high- spin and low- 
spin fep cluster at P~ 120 GPa. Gontributions of different 
vibrational modes to (3 are given by impulses. 



ter, this mode occurs at somewhat lower frequencies, 
and is spread out among several similar modes. This 
mode does not fractionate iron as intensely in the high- 
spin cluster as it does in the low-spin cluster; P is nearly 
four per mil lower than the low-spin cluster after pass- 
ing through these vibrational modes. The mode at 560.3 
(559.8) cm~^ in the low-spin cluster, mostly involving 
motion of the equatorial Fe04 plane against the oxygens 
at the apex, also makes a much larger contribution to P 
than the corresponding vibration in the high-spin clus- 
ter. The next two vibrational modes make roughly equal 
contributions to P in both the high- and low-spin clus- 
ters. The last two modes, involving primary stretching 
of the FeOe octahedron, with the iron in phase with the 
magnesium atoms, again make stronger contributions for 
the low-spin cluster than for the high-spin cluster. 

The magnitudes of the increases in p with pressure 
has implications for the isotopic signatures of accreting 
planetary bodies. First, an individual planet in isotopic 
equilibrium should be isotopically stratified with heavier 
isotopes progressively concentrated at depth within the 
planet. Such a stratification could possibly be established 
during a "magma ocean" stage as iron isotopes sample 
coordination environments in liquid silicates that would 
be similar to fep or fepv. An interesting question is 
whether a spin transition would take place in the melt. 
If coordination environments are structurally closer to 
fepv [2F, the spin transition may be supressed and the 
isotopic signature could be lighter than a solid of the 
same composition. Freezing of a magma ocean might be 
expected to give rise to an isotopically stratified mantle 
with the heavy isotopes accumulating preferentially at 
the base, depending on the relative densities of the liquid 
and solid phases. 

If an equilibrium isotopic profile, with ^^Fe/^^Fe in- 
creasing with depth, can be established after impact 
events, evaporative loss of surface material during ex- 
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FIG. 3: Major fractionating vibrational modes shown in 
Fig. 2. Numbers in italics give vibrational frequencies in 
cm~^. The modes at 560 cm~^ and 593 cm~^ involve similar 
displacements of the central iron atom. 



treme heating would be isotopically light, thus the larger 
planets should be isotopically heavier than the smaller 
planets, because the larger planets can shield greater re- 
serviors of material at high-pressure where the /3 fac- 
tor is larger. This may be part of the reason why the 
Earth-Moon system appears to be isotopically heavier 
than Mars and Vesta [26 . 

The pressure dependence of P is sensitive to crystal- 
lographic environment; /3fepv increases more than a fac- 
tor of two over high-spin /3fep, probably because of the 
greater compressibility of the coordination shell in fepv. 
However, the increase is compensated completely at the 
core-mantle boundary by the spin transition in fep. The 
similarity between /3fepv and /3fep indicates little differen- 
tial isotopic fractionation between these phases after the 
spin-transition takes place in the lower mantle. Above 
this transition, we expect fepv to be enriched in ^^Fe 
over fep by a factor of two. 
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